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ABSTRACT 

Recent Chandra observations of an outflowing gas in GRO J 1655^0 resulted in a suggestion by Miller et 
al. (2006) that the wind in this system must be powered by a magnetic process that can also drive accretion 
through the disk around the black hole. The alternative explanations, of radiation pressure or thermally driven 
flows, were considered unsatisfactory because of the highly ionized level of the gas and because of the derived 
small distance from the black hole, well inside the minimum distance required for an efficient X-ray heated 
wind. The present paper shows that there is a simple photoionized wind solution for this system where the gas 
is much further out than assumed by Miller et al., at r/rg — 10"*'^^^^. The expected wind velocity, as well as 
the computed equivalent widths of more than 50 absorption lines in this single-component ID model, are all in 
good agreement with the Chandra observations. 

Subject headings: binaries: accretion, accretion disks - spectroscopic - black hole physics - X-rays: stars - 
stars: individual (GRO J1655-40) 



1. INTRODUCTION 

High resolution spectroscopy of transient X-ray sources is 
a powerful tool for studying the nature of such systems. The 
recent paper by Miller et al. (2006, hereafter M06) is an ex- 
cellent example that illustrates the power of such methods. 
These authors analyzed the high energy transmission grating 
■ (HETG) Chandra observations of GRO J 1655-40 (Orosz & 
Baily 1997), a stellar-mass binary system containing a 7 Mq 
BH and a F3IV-F6IV companion with a mass of 2.3 Mq at a 
mean distance of about 10^ gravitational radii {rg = GM/c^). 
The observations were obtained during an X-ray bright phase 
when the intrinsic luminosity of the source approached 3.3 x 
10^^ erg s^', equivalent to L/LEdd^ 0.03. They show the 
clear signature of an outflow ("X-ray wind"). M06 measured 
the equivalent widths (EWs) and the line widths of more than 
70 identifiable X-ray lines. Using empirical fits to the line 
profile, they derived column densities for the more abundant 
ions (mostly H-like and He-like lines of oxygen, neon, mag- 
nesium, silicon, sulfur, argon, calcium, iron and nickel) and 
obtained reliable estimates of the gas composition and level of 
ionization. Fitting the unblended lines indicate outflow veloc- 
ities of 300-1000 km s"' and line widths of 300-500 km s"'. 
This data set is not the first of its kind (e.g. Miller et al. 2004) 
but, so far, it is the highest quality one. 

According to M06, the X-ray wind in GRO J 1655-40 must 
be powered by magnetic processes (pressure generated by 
magnetic viscosity or magnetocentrifugal forces, e.g. Bland- 
ford and Payne 1982). The reason for favoring this interpre- 
tation is the inefficiency of two alternative mechanisms, ra- 
diation pressure force and thermal expansion, to explain the 
observed properties. This conclusion, claimed to be the best 
evidence so far for a magnetically driven disk wind, was based 
on a detailed photoionization model for the X-ray heated gas. 
The outflowing gas in this model is very close to the BH, 
some 450ri, from it, where the escape velocity is about 2 x 10^ 
km s^'. This location was inferred from trying to match the 
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observed EWs of more than 70 absorption lines. The best fit 
model suggests a very high ionization parameter, gas density 



,24.15 



of n ~ 5.6 X 10 cm^ and a column density of Nh — 10 
cm^^. The observed wind velocity is more than an order of 
magnitude smaller than the escape velocity at the suggested 
location. This was explained to be due to the high inclination 
of the disk (67-85 degrees) in this system. The highest veloc- 
ity component is perpendicular to the accretion disk and the 
observed velocities are much smaller. Given this location and 
level of ionization, the radiation pressure force is too small 
to drive the wind. A thermal wind solution was rejected on 
grounds that the gas location is orders of magnitude inside the 
region where such winds are likely to develop (e.g. Begelman, 
McKee and Shields 1983). The rejection of these two mecha- 
nisms leave magnetic processes as the only viable explanation 
for the outflow. 

In this paper I show that there is, in fact, a simple wind 
solution that is consistent with all the observed properties of 
this system. §2 presents this solution and §3 discusses several 
other implications. 

2. A THERMAL WIND MODEL FOR GRO 11655^0 
2.1. General considerations 

Thermal winds are expected to arise in various astro- 
physical situations involving X-ray irradiated accretion disks. 
These includes disks around galactic X-ray sources as well as 
the inner regions of active galactic nuclei (AGNs). In such 
cases, the central X-ray source provides the energy to lift the 
gas from the surface of the disk at large distances and to heat 
it to a high enough temperature that enables escape via an ex- 
panding wind. Such winds have been discussed extensively 
in the literature. Here I focus on the Begelman et al. (1983) 
analytical approach and the follow-up 2D time-dependent nu- 
merical calculations of Woods et al. (1996, hereafter W96). 

Important parameters in X-ray irradiated thermal wind 
models are the Compton temperature of the gas, 7^, and the 
Compton radius, r^, which is the distance at which the escape 
velocity equals the isothermal sound speed at the Compton 
temperature. This radius is given by 

9.8 X \(fM{BH) 
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where Tcs = 7c/10^K. According to the theoretical esti- 
mates, X-ray heated winds can develop at any r > 0.1 rp. The 
more accurate numerical calculations (e.g. Fig. 4 in W96) 
give even smaller radii, as small as O.Olrc. There are var- 
ious wind regions depending on the luminosity (isothermal 
wind, steadily heated free wind and gravity inhibited wind, 
see W96 for more details). The Compton temperature of the 
gas depends on the spectral energy distribution (SED) of the 
continuum source. In the case of GRO 11655^0, this is 
obtained directly from the observations, after correcting for 
galactic absorption. There is little uncertainty in the chosen 
SED since most of the energy is emitted over the observed 
Chandra band. Here I adopt the SED described in M06 which 
is a combination of a thermal source with ^7 = 1.35 keV and a 
high energy power-law of the form Le °^ E^'^-^'^. This results 
in 7c ~ 1.4 X 10^ K. Thus, rc — lO'^'^ cm and an effective 
thermal wind is likely to develop at all distances exceeding 
about 10^*^^ cm (lO^^-^r^) and perhaps much smaller. 

The simple ID model presented here assumes a point-like 
X-ray source in the center of the disk and a An ionized gas dis- 
tribution that may have different density and column density 
at different polar angles. This is not a full numerical model but 
rather a demonstration that a differentially expanding, pho- 
toionized thermal wind can explain the Chandra observations 
of GRO 11655^0. Reahstic wind models combine the micro- 
physics of the ionized gas with the gas dynamics, in a time de- 
pendent fashion. Such calculations are are beyond the scope 
of the present work or the work of M06. However, as shown 
by Begelman et al. (1983), many of the wind properties can 
be inferred from general considerations. Further confirmation 
can be found in Chelouche and Netzer (2005) who describe a 
phenomenological approach of a polytropic wind that shares 
many properties with numerical ID and 2D wind simulations. 
Simple estimates of this type are nevertheless limited in the 
abUity to estimate the wind velocity and density at various 
locations. 

Thermal wind are characterized by two regimes, before and 
after crossing the critical (sonic) point. The present work ad- 
dresses the region around and beyond the sonic point, where 
most of the observed line opacity is assumed to take place. In 
this region the velocity increases slowly with radius in a way 
that can be expressed as 

Vwind °^r^ , (2) 

where in most cases of interest < x < 0.5. The wind speed at 
r^rc approaches 'i-A time the sound speed in the fully ion- 
ized gas. AS shown by W96, thermal outflows can develop in 
gas whose temperature is much below Tc (e.g. W96 Fig. 7). 
Such gas is carried away with the flow and can reach, under 
some circumstances, its Compton temperature. The case ad- 
dressed here is that of the cooler gas although, for simplicity, 
some of the estimates are based on the assumption T^ylnd = Tc- 
Continuity in a spherical geometry (see however comments 
in §3) requires that 

„ oc . (3) 

Assign no as the gas (hydrogen) density at the base of the 
wind, just beyond the sonic point, where the distance from 
the BH is tq. This gives a maximum column density of 

NH=noro/il+x). (4) 

The ionization parameter of the gas is given by 

^ = ^, (5) 



where Lx is the integrated luminosity above 1 Rydberg and 
where at the base of the wind = Lx/nor^. Thus, <^ cx r^. 

Experimenting with the above SED and a detailed pho- 
toionization code (see §2.2), and changing the ionization pa- 
rameter until it results in a fractional ionization that fits the 
observations to within a factor two, I obtain ^ ~ lO''. This is 
15-50 times smaller than the value suggested by M06^. Given 
the column densities inferred from the EWs of several unsat- 
urated Unes (M06 supplementary Table 1), and assuming near 
solar composition, I derive a minimum required column den- 
sity of Nh — 10^^'^ cm^^, about a factor 3 smaller than the 
column density used in M06. The combination of and Nh, 
together with Lx, can be used to solve for no and the max- 
imum value of ro which is about lO'^ cm (~ lO^r^). This 
location is outside the minimum distance (O.lrc or smaller) 
required for a thermal wind. Some 90% of the wind opac- 
ity is achieved over the distance 1-lOro. For GRO 11655^0, 
this corresponds to 10"*^ ^ r/rg ^ 10^-^. These dimensions 
fit nicely with the expected dimensions of the accretion disk 
in this system. The expected wind velocities over this region 
are similar to the velocities inferred from the observed line 
blueshifts in the spectrum of the source. 

Thus, estimates based on the observed ionization, column 
density and velocity of the gas in GRO 11655^0, all agree 
with the conditions expected in a thermal wind around an X- 
ray irradiated accretion disk. 

2.2. A photoionized wind model for GRO 31655^0 

Detailed photoionization calculations were carried out to 
compare the predicted spectra, under such conditions, with 
the Chandra observations. The photoionization code used, 
ION06, is the 2006 version of the code ION most recently 
described in Netzer (2004) and in Netzer et al. (2005). The 
only significant modification is the introduction of new radia- 
tive and di-electronic recombination rates described in Bad- 
nell (2006) and references therein. The latest atomic data for 
iron L-shell lines, published by Landi and Gu (2006), were 
used in order to faciUtate a meaningful comparison with the 
large number of Fe xxii and Fe xxiii lines observed in the 
spectrum of this source. The SED and Lx are those given in 
M06. The gas composition was assumed to be solar except 
for calcium and iron where, following the suggestion of M06, 
twice the solar abundances were assumed. 

Several velocity profiles were tested assuming < x < 0.5. 
Line widths of 250-500 km s^' were chosen to agree with 
the observations. These velocities can be interpret as due to 
internal (sound or turbulent) or global (wind expansion) mo- 
tion. The results are not sensitive to the exact values of t,, 
ro and no provided they are within a factor ~ 2 of the values 
given above. Several models with ^ decreasing with r have 
also been computed and found to be of similar quality. 

Two additional parameters require further explanation. The 
first is the covering fraction of the outfiowing gas which, given 
the disk-wind geometry, may not be the same for the emis- 
sion and absorption lines. As noted by M06, the very weak 
emission lines and the lack of clear P-Cygni profiles suggest 
a small emission measure of the X-ray gas. This can be inter- 
preted as different density and/or colunm density at different 
polar angles. For example, the W96 numerical calculations 
(Fig. 7 and Tables 2-4) clearly show that column densities 
along the plane of the disk can be significantly larger than the 

^ The reason for this discrepancy, which is a major cause for the very 
different model presented here, is not clear 
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TABLE 1 

Parameters of thermal wind models. 



^0 rp 

1.1 X 10^ 10'" " cm 10" " 10^3.73 ^^-2 q 3qq 0.99±0.40 

9.2 X 10^ lO'"-*' cm 10" 25 cm-^ 10^3 cm'^ 0.3 400 1.03 ±0.46 



column perpendicular to the disk surface. Such a scenario will 
produce weak emission lines. In the ID photoionization cal- 
culations presented here, this is taken into account by assum- 
ing a small emission (global) covering factor, Q.^m/'^'^ < 0.2. 

There are indications in the spectrum that the absorption 
(line-of-sight) covering factor is also smaller than unity. This 
can be the result of a patchy absorber or due to scattering of 
the X-ray continuum into the Une of sight. For example, sev- 
eral of the longer wavelength lines, e.g. those of neon and 
magnesium, are resolved yet the optical depths inferred from 
the profiles of the low and the high series lines are inconsistent 
and suggest an unobscured continuum. A second example is 
the pair of Fe xxiv lines at 10.67,10.62A. The optical depths 
of the two differ by a factor 2 yet the observations suggest al- 
most identical depths. This indicates saturation in both lines 
combined with an unobscured continuum. A rough estimate 
based on visual inspection suggests an absorption covering 
factor of about 0.75 . The EWs Usted in M06 must be corrected 
for this factor in order to compare with the calculations. 

Regarding the gas density, in principal this can be di- 
rectly obtained from the relative intensity of the Fe Xll lines 
since one of the fines, at 11.92A, originates from an excited 
2s^2p-^Py2 level at 14.6 eV above ground. The known col- 
lision strength and decay rate of this level indicate a critical 
density of order lO'^ cm^^ for the temperature found here 
(~ 2 X 10^ K). However, this is only true for optically thin 
transitions while for the density estimated here (about 10^^ 
cm^^ at the base of the wind) and the assumed column den- 
sity, the lines in question are optically thick (the ground level 
transitions are starting to saturate while the optical depth of 
the II.92A fine is of order unity). This changes the line ra- 
tios in favor of the smallest optical depth line and results in a 
satisfactory agreement between the observations and the cal- 
culations of the EWs of all Fe XXII lines. 

The results of two models with different assumptions on the 
radial velocity of the wind, x = and x = 0.3, are summarized 
in Table 1 . The calculated spectrum for the x = atmosphere 
is shown in Fig. 1 on a scale similar to the one shown in M06 
(supplementary Fig. 1). The computed spectrum was con- 
volved with the HEG instrumental profile to enable a direct 
comparison with the observations. The diagram shows the 
great similarity of the model with the Chandra observations'*. 

To make a more quantitative comparison I show, in Fig. 
2, a comparison of 53 EWs of observed and calculated satu- 
rated and unsaturated fines. These lines are basically the ones 
used by M06 in their Fig. 2 except that the present calcu- 

^ Note that the flux scale used here corresponds to the absorption corrected 
spectrum while M06 plot the observed spectrum. Note also that the Fe XXIV 
L-edge at 6.06A looks artificially prominent because the current model does 
not include the high order L-shell lines of this ion and the Ni XXVI line at 
6.05A 
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Fig. 1 . — A thermal wind model for GRO 11655^0 (the absorption spec- 
tram of model A). The continuum is corrected for galactic absorption and an 
absorption covering fraction of 0.75 is assumed. 



lations do not include nickel and sodium lines. The values 
plotted are EW(observed)/EW(model) where EW(model) is 
obtained by assuming Gaussian fine profiles (damping winds 
are not important at such velocities) and where the computed 
EWs are corrected for the blending of several line pairs. The 
diagram shows a uniform distribution with no obvious ion- 
ization trends. The mean and the standard deviation of this 
ratio, for the two models detailed in Table 1, are 0.99 ± 0.40 
for the X = model and 1 .03 ± 0.46 for the x = 0.3 model. 
Both are considered satisfactory given the measurement un- 
certainties and the simpficity of the model and both are of 
comparable quality to the M06 fit (their Fig. 2). There are 
only two notable deviations. The first is Ne IX 13.45A where 
EW(observed)/EW(model)=2.3±0.5. This is suspected to be 
due to an overestimate of the measured EW and can be tested 
by measuring a different Ne ix line observed in the spec- 
trum but not listed in M06. The second is Fe xxii 11.77A 
where EW(observed)/EW(model)=0.43±0.06. However, the 
spectrum shows several additional Fe XXII lines that were not 
identified by M06 (e.g. at 11.54 and 11.42A) that seem to 
be in better agreement with the model, although their EWs 
was only estimated from the plot and hence was not included 
in the comparison. The ground-to-excited level line ratios in 
Fe XXll is in good agreement with the observations. 

3. DISCUSSION 

The Chandra observations of GRO J 1655-40 are perhaps 
the best indication, so far, of X-ray heated flows in galactic 
sources containing accretion disks. M06 have argued that the 
origin of this wind is magnetic because two alternative mech- 
anisms that can drive such flows are inefficient. One possibil- 
ity, acceleration by radiation pressure force, was ruled out on 
ground that the force multipfier in such highly ioruzed gas is 
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Fig. 2. — Observed over calculated EWs for the absorption lines in model 
A. FviU squares represent lines from H-like ions, empty squares are lines of 
He-like ions and triangles are lines of Fe xxil, Fe xxill and Fe xxiv. 



too small to drive the observed flow. The models presented 
here agree with this conclusion. The value of the force mul- 
tiplier (~ 1.5) is more than an order of magnitude too small 
to drive a wind given L/LEdd of the source. The second pos- 
sibiUty, of an X-ray heated wind, was rejected on the basis of 
the derived location of the absorbing gas which was obtained 
by computing a detailed photoionization model to fit the ob- 
served spectrum. This was found by M06 to be well inside the 
minimum distance where thermally driven winds are thought 
to be important. 

The model presented here is of similar quality to the M06 
model yet the gas location is some two orders of magnitude 
further out. This is the region where thermal. X-ray heated 
outflow are likely to be important. While not a complete 2D 
time dependent wind model, all the properties found here, in- 
cluding the location, velocity and column density of the gas, 
are consistent with theoretical expectations. The mass loss 
rate in this model is m„ ~ 2 x 10'^i2gm/4;r gr s^' (note a 
comment below about non-radial streamlines). For an accre- 
tion efficiency of 0. 1 , this is about a factor of 5QQ.em/ 47r larger 



than the accretion rate by the central object. Given the small 
derived £lem, this is in general agreement with the Begelman 
etal. (1983) results. 

Realistic 2D numerical simulations that include the possi- 
biUty of a non-uniform medium and the disk- wind interaction 
are required to fully explain the spectrum of GRO J 1655-40. 
Such calculations must address several important issues that 
could not be included in the present model: 1 . Thermal insta- 
bilities in the outflowing gas. 2. A multi-component patchy 
absorber. 3. Non-radial streamlines that will modify the sim- 
plified continuity condition (eqn. 3) and hence also the radial 
dependence of (e.g. a decrease of ^ with r). 4. Radial 
dependence of the covering factor. 5. Time dependent Lx 
reflecting the observed variability of the source. The simple 
ID model presented here is already in good agreement with 
the observations and is likely to be a major component of the 
more advanced computations. 

All the above does not exclude the possibility that magnet- 
ically driven winds are important sources of driving material 
near galactic and extragalactic BHs, including GRO J1655- 
40. Such winds have been proposed and explored in numerous 
papers and provide a natural scenario for many observed sys- 
tems (e.g. Blandford & Payne 1982; Hawley, Gammie & Bal- 
bus, 1995). However, given the success of the simple thermal 
wind model to explain the observations of GRO 11655^0, 
there seems to be no need to invoke a magnetically powered 
wind as a major contributer in this case. 
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